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We address the long standing problem of the dependence of the electrostatic persistence length 
l e of a flexible polyelectrolyte (PE) on the screening length r s of the solution within the linear 
Debye-Hiickel theory. The standard Odijk, Skolnick and Fixman (OSF) theory suggests l e oc rf, 
while some variational theories and computer simulations suggest l e oc r s . In this paper, we use 
Monte-Carlo simulations to study the conformation of a simple polyelectrolyte. Using four times 
longer PEs than in previous simulations and refined methods for the treatment of the simulation 
data, we show that the results are consistent with the OSF dependence l e oc r\. The linear charge 
density of the PE which enters in the coefficient of this dependence is properly renormalized to take 
into account local fluctuations. 
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I. INTRODUCTION 

Despite numerous theoretical studies of polyelectrolyte 
(PE), due to the long range nature of the Coulomb in- 
teraction, the description of their conformation is still 
not as satisfactory as that of neutral polymers. One 
of the longest standing problem is related to the elec- 
trostatic effect on the rigidity of a PE. In a water so- 
lution with monovalent ions, within the Debye-Hiickel 
linear screening theory, the electrostatic interaction be- 
tween PE charged monomers has the form: 
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where r is the distance between monomers, D is the di- 
electric constant of water, e is the elementary charge, and 
r s is the Debye-Hiickel screening length, which is related 
to the ionic strength / of the solution by r 2 = AttIbI- 
(Ib = e 2 /DksT is the Bjerrum length, T is the temper- 
ature of the solution). 

The rigidity of a polymer is usually characterized by 
one parameter, the so called persistence length l p . For 
a polyelectrolyte chain, besides the intrinsic persistence 
length Iq which results from the specific chemical struc- 
ture of the monomers and bonds between them, the total 
persistence length also includes an "electrostatic" contri- 
bution l e which results from the screened Coulomb inter- 
actions between monomers: 



lp — lo + h ■ 



(2) 



Because the interaction (]l]) is exponentially screened 
at distances larger than r s , early works concerning the 
structure of the PE assumed that l e is of the order of r s . 
However, this simple assumption was challenged by tha 
pioneering works of Odijkcl and Skolnick and Fixmancl 
(OSF), who showed that Debye-Hiickel interaction can 
induce a rod-like conformation at length scales much 



larger than r s . Their calculation gives 
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where rjo is the linear charge density of the PE. Because 
l e oc r 2 s , it can be much larger than r s at weak screening 
(large r s ). 

Although the idea that electrostatic interaction en- 
hances the stiffness of a PE is qualitatively accepted and 
confirmed in many experiments, the quadratic depen- 
dence of l e on the screening length r s is still the subject 
of many discussions. In the work of OSF, the bond angle 
deflection was assumed to be small everywhere along the 
chain, what is valid for large lo. They suggested that if lo 
is not small but l e is large enough (week screening) , their 
assumption is still valid. Ref. 0, however, has questioned 
this assumption especially when lo is so small that the 
bond angle deflection is large before electrostatics comes 
into play and rigidities the chain. 

A significant progress was made by Khokhlov and 
Khachaturian (KK) who proposed a generalized OSF 
theoryQ for the case of flexible polyelectrolyte (small ^o)- 
It is known that in the absence of screening (r s — > oo), 
the structure of a polyelectrolyte can be conveniently de- 
scribed by introducing the concept of electrostatic blobs. 
A blob is a chain subunit within which the electrostatic 
interaction is only a weak perturbation. The blob size 
£ is related to the number of Kuhn segments g within 
one blob as £ = log 1 ! 2 ■ The condition of weak Coulomb 
interaction suggest that the electrostatic self energy of 
a blob, (rjoglo) 2 1 is of the order of ksT. This leads 
to £ ~ (DkBTlo/rjo) 1 ^ 3 - At length scale greater than £, 
Coulomb interaction plays important role and the string 
of blobs assumes a rod-like conformation, with the end- 
to-end distance proportional to the number of blobs. 

Using this blob picture, KK proposed that OSF theory 
is still applicable for a flexible PE provided one deals 
with the chain of blobs instead of the original chain of 
monomers. This means, in Eq. (^|), one replaces the 
bare linear charge density 770 by that of the blob chain 
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i] = rjoglo/t;. The intrinsic persistence length Iq should 
also be replaced by £. As a result, the total persistence 
length of the flexible PE reads: 
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Thus, in KK theory, despite the flexibility of the PE, its 
electrostatic persistence length remains quadratic in r s . 
Small Iq only renormalizes the linear charge density from 
Vo to V- 

Note that r s is implicitly assumed to be larger than the 
blob size £ in KK theory (weak screening). For strong 
screening r s < £, there are no electrostatic rigidity and 
the chain behaves as flexible chain with the Debye-Huckcl 
short range interaction playing the role of an additional 
excluded volume interaction. 

A number of variational calculations have also been 
proposed to describe more quantitatively the structure 
of flexible chain. These calculations, although based on 
different ansatz, have the same basic idea of describing 
the flexible charged chain by some model of noninteract- 
ing semiflexible chain and variationally optimizing the 
persistence length of the noninteracting system. Surpris- 
ingly, while some of these calculations support the OSF- 



KK dependence l e oc r 2 such as Refs. ^gflj, other c. 



lations found that l e scales linearly with r s insteac 
However, because variational calculation results depend 
strongly on the variational model Hamiltonian, none of 
these results can be considfpsd-xoncjjusive. 

Computer simulations^OOEjotj also have been 
used to determine the dependence l e on r s and to ver- 
ify OSF or variational theories. Some of these papers 
claim to support the linear dependence of l p on r s . The 
simulation of Ref. [l5| concludes that the dependence of l p 
on r s is sublinear. Thus, the problem of the dependence 
l e (r s ), despite being very clearly stated, still remains un- 
solved for a flexible PE. More details about the present 
status of this problem can be found in Ref. [r| 

In this paper, we again use computer simulations to 
study the dependence of l e on r s . The longest poly- 
electrolyte simulated in our paper contains 4096 charged 
monomers, four times more than those studied in previ- 
ous simulations. This allows for better studying of size ef- 
fect on the simulation result. Furthermore, we use a more 
refined analysis of the simulation result, which takes into 
account local fluctuations in the chain at short distance 
scale. Our results show that OSF formula quantitatively 
describes the structure of a polyelectrolyte. 

The paper is organized as follows. The procedure of 
Monte-Carlo simulation of a polyelectrolyte using the 
primitive freely jointed beads is described in the next sec- 
tion. The data for the end-to-end distance R ee is given. 
In Sec. Ill, we analyze this data using the scaling argu- 
ment to show that it is consistent with OSF theory. In 
Sec. IV, we analyze the data for the case of large r s , 
where excluded volume effect is not important, in order 
to extract l e and again show that it obeys OSF theory 
in this limit. In Sec. V, we use the bond angle correla- 



tion function to calculate l e and to confirm the result of 
Sec. IV. The good agreement between l e calculated us- 
ing different methods further suggests that OSF theory 
is correct in describing a polyelectrolyte structure. We 
conclude in Sec. |Vl| . 

Several days after the submission of our papeii-to the 
Los Alamos preprint archivetll, another paperEj with 
Monte-Carlo simulations for PE molecules in the same 
range of lengths appears in the same archive. Results of 
this paper are in good agreement with our Sec. III. 



II. MONTE-CARLO SIMULATION 

The polyelectrolyte is modeled as a chain of N freely 
jointed hard spherical beads each with charge e. The 
bond length of the PE is fixed and equal to Ib, where 
Ib = e 2 / DksT is the Bjerrum length which is about 
7A at room temperature in water solution. Thus the 
bare linear charge density of our polyelectrolyte is 770 = 
e/ls- Because we are concerned about the electrostatic 
persistence length only, the bead radius is set to zero so 
that all excluded volume of monomers is provided by the 
screened Coulomb interaction between them only. For 
convenience, the middle bead is fixed in space. 

To relax, the PE configuration globally, the pivot 
algorithmic is used. In this algorithm, in an attempted 
move, a part of the chain from a randomly chosen 
monomer to one end of the chain is rotated by a ran- 
dom angle about a random axis. This algorithm is known 
to be very efficient. A new independent sample can be 
produced in a computer time of the order of N , or in 
other words, uncorrelated samples are obtained every few 
Monte Carlo (MC) steps (one MC step is defined as the 
number of elementary moves such that, on average, ev- 
ery particle attempts to move once). To relax the PE 
configuration locally, the flip algorithm is used. In this 
algorithm, a randomly chosen monomer is rotated by a 
random angle about the axis connecting its two neigh- 
bor (if it is one of the end monomers, its new position is 
chosen randomly on the surface of a sphere with radius 
Ib centered at its neighbor.) In a simulation, the num- 
ber of pivot moves is about 30% of the total number of 
moves. The usual Metropolis algorithm is used to accept 
or reject the move. About 1 -r 2 x 10 4 MC steps are run 
for each set of parameters (N, r s ), of which 512 initial 
MC steps are discarded and the rest is used for statistical 
average (due to time constrain, for TV = 4096, only 2000 
MC steps are used). Two different initial configurations, 
a Gaussian coil and a straight rod, were used to ensure 
that final states are indistinguishable and the systems 
reaches equilibrium. 

The simulation result for the end-to-end distance R ee 
of a polyelectrolyte for different N is plotted in Fig. |l| as 
a function of the screening radius r s of the solution. At 
very small r s , Coulomb interactions between monomers 
are strongly screened and the chain behaves as a neutral 
Gaussian chain with R ee = IbVN — 1. At very large 
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FIG. 1: The square of the end-to-end distance of a polyelec- 
trolyte R^ e as a function of the screening length r s for chains 
with different number of monomers N: 64(o), 128(+), 256 
(□), 512 (x), 1024 (A), 2048(*), and 4096 (♦). The arrows on 
the right side show B^, e obtained using unscreened Coulomb 
potential V(r) = e/r. 



r s 3> N, Coulomb interactions between the monomers 
are not screened and R ee is saturated and equal to that 
of an unscreened PE with the same number of monomers 
(see the arrows in Fig. |l|). 

Three different methods are used to verify the validity 
of OSF theory for flexible PE: i) study of the scaling de- 
pendence of R ee on r s in whole range of r s , ii) extraction 
of l e in the large r s limit and iii) analysis of the bond cor- 
relation function. In the next three sections, we discuss 
these methods in details together with their limitations. 
Comparison with previous simulations is also made to ex- 
plain their results which so far have not supported cither 
of the theories. 



III. SCALING DEPENDENCE OF R ee ON r s 

Let us first describe theoretically how the chain size 
should behave as a function of the screening radius r s 
when r s increases from to oo. 

When r s <C Ib, the Coulomb interaction is strongly 
screened. Because there are no other interaction present 
in our chain model of freely jointed beads, the chain 
statistic is Gaussian. Its end-to-end distance R ee is pro- 
portional to the the square root of the number of bonds 
and independent on r s : 
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When r s ^> Is, the chain persistence length is domi- 
nated by the Coulomb contribution l p ~ Z e . If N is very 
large such that the chain contour length NIb is much 
larger than l e then the chain behaves as a linear chain 



with Nls/le segments of length l e each and thickness r s . 
The excluded volume between segments is v ~ l 2 r s , and 
the end-to-end distanced: 
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At larger r s where l e becomes comparable to the PE 
contour length, the excluded volume effect is not impor- 
tant. In this case, the chain statistics is again Gaussian 
and 
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Finally, at even larger r s when l p is greater than Nl b , 
the chain becomes a straight rod with length independent 
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If l e oc rf, the transition from the scaling range of 
Eq. (§) to Eq. (@) happens at r s ~ IbN 1 / 4 , while the 
transition from the scaling range of Eq. (^) to Eq. (||) 
happens at r s ~ IbN 1 / 2 . On the other hand, if l e oc r s , 
both transitions from the scaling range of Eq. (||) to Eq. 
(0) and from the scaling range of Eq. (0) to Eq. (g) 
happen at r s ~ IbN . This means, there is no scaling 
range of Eq. (0) in this theory. 

Thus, one can distinguish between the OSF result, l p oc 
r 2 s , and the variational result, l p oc r s by plotting the 
exponent a = 91n[i?g e ]/91nr s as a function of lnr s . The 
schematic figure of this plot is shown in Fig. ||. OSF 
theory gives plateaus at a = 6/5 and 2, and when r s > 
IbN 1 / 2 , a drops back to 0. Variational theories, on the 
other hand, would suggest one large plateau at a = 4/5 
up to r s ~ IbN. 

The simulation results for a are shown in Fig. [|for dif- 
ferent TV. One can see that as TV increases, the agreement 
with OSF theory becomes more visible. Note that the 
plateaus in Fig. Bare scaling ranges, and relatively sharp 
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FIG. 3: Simulation result for a as a function of r s for different 
N: 64(o), 128(+), 256 (□), 512 (x), 1024 (A), 2048(*), and 
4096 (♦). They agree reasonably well with the solid curve of 
Fig. 0, suggesting that OSF theory is correct. 



transitions between plateaus are valid only for TV — > oo. 
For a finite TV, the plateaus may be too narrow to be ob- 
served and can be masked in the transition regions. This 
explains why one cannot see the plateau at a = 2 in our 
results. Nevertheless, the tendencies of a to develop a 
plateau at a = 6/5, then to grow higher toward a = 2 at 
larger r s and finally to collapse to zero when approach- 
ing relatively small r s — IbV~N are clearly seen for large 
TV. Thus, generally speaking, the curves agree with OSF 
theory much better than with variational theories (where 
a is supposed to be about 4/5 and to decrease to zero 
only when r s — > IbN , i.e. at much larger r s than what 
observed) . 

Fig. ^ also shows one reason why similar simulations 
done by other groups do not support OSF theory. All of 
these simulations are limited to 512 charges. As one can 
see from Fig. [| the curves for TV < 512 do not permit 
to discriminate between the two theories as clearly as the 
case TV — 2048 or 4096. Only when TV becomes very large 
can scaling ranges with a > 1 show up and one observes 
better agreement with OSF result. 



IV. LARGE 



LIMIT 



In this section, we attempt to extract directly from 
the simulation data the persistence length in order to 
compare with OSF theory. To do this, one notices that 
even a chain with excluded volume interaction behaves 
as a Gaussian chain when its contour length is very short 
such that it contains only a few Kuhn segmented In this 
case, one can use the Bresler-Frenkel formulacj to de- 
scribe the relationship between the end-to-end distance 



R ee and the chain persistence length l p : 

R 2 ee = 2Ll p - 2l 2 p [l - exp(-L/i p )] 



(9) 



where L is the contour length of the chain. 

For our polyelectrolyte, this formula can be used for 
large r s when the persistence length is of the order of 
R ee or larger. However, one cannot use the bare contour 
length Lq = (TV — 1)1 b in the Eq. (|J) because the chain 
where OSF theory is supposed to be applicable is not the 
bare chain but an effective chain which takes into account 
local fluctuations. The contour length L of this effective 
chain is 



L = Ne/rj 



(10) 



where r\ is the renormalized linear charge density of the 
PE. 

In KK theory, the effective chain is the chain of elec- 
trostatic blobs, and the normalized charge density is 
V = Voglo/S,- However, the standard blob picture can 
only be used to describe flexible weakly charged chains 
where the fraction of charged monomers is small so that 
the number of monomers, <?, within one blob is large 
and Gaussian statistics can be used to relate its size 
and molecular weight. Because, for a given number of 
charged monomers, Monte-Carlo simulation for weakly 
charged polyelectrolyte is extremely time consuming, all 
monomers of our simulated polyelectrolyte are charged. 
In this case, the neighbor-neighbor monomers interaction 
equals ksT. This makes g — 1 and the standard picture 
of Gaussian blobs does not apply. Thus, in order to treat 
our data, we assume that both l p and r\ are unknown 
quantities. 

To proceed further, one needs an equation relating r\ 
and l p , and in order to verify OSF theory, we could use 
their formula 



l p = r?rl/4Dk B T , 



(11) 



for this purpose. Thus, we could substitute Eq. jlO| ) and 
( pd| ) into Eq. (^), and solve for rj using R 2 e obtained 
from simulation. If OSF theory is valid, the obtained 
values of r\ should be a very slow changing function of 
r s . In addition, in the limit TV — > oo, they should also be 
independent on TV. 

The OSF equation (|J), however, was derived for the 
case r s -C L while in our simulation, the ratio r s /L is 
not always small. Therefore, instead of Eq. n ([Tl|), we use 
the more general Odijk's finite size formulaEI 
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for the persistence length l p . When L » r Sl the term in 
the square brackets is equal to 3 and the standard OSF 
result is recovered. On the other hand, when r s L, 
the persistence length l p saturates at rj 2 L 2 /72DkBT. 

Below, we treat our Monte-Carlo simulation data with 
the help of Eq. @ using Eq. Q and @ for L and 
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FIG. 4: The linear charge density r\ as a function of the screen- 
ing length for different AT: 64(o), 128(+), 256 (□), 512 (x), 
1024 (A) and 2048(*). The thick solid line is the theoretical 
estimate which is the numerical solution to Eq. ([l3|) , ( |l4| ) and 
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FIG. 5: Plots of L/r s as a function of r s calculated with the 
help of Eq. (§), @ and using our data for ( + ) an d 
using unperturbed r\ — T]o as in Ref. [TB] (o). The chain with 
N = 1024 is used. 



l p . The results for rj are plotted in Fig. || for different 
PE sizes N. As one can see, at large r s , r\ changes very 
slowly with r s , and as N increases, tends to saturate at 
an N independent value. 

It should be noted that the lines r)(r s ) in Fig. [| un- 
physically start to drop below certain values of r s . This is 
because at smaller r 8 , the electrostatics-induced excluded 
volume interactions between monomers become so strong 
that the right hand side of Eq. (||) (which is derived for a 
Gaussian worm like chain) strongly underestimates R ee . 

Even though the picture of Gaussian blobs does not 
work for our chain, rj can still be calculated analytically 
in the limit L r s (N — > oo). Indeed, let us assume 
that the effective chain is straight at length scale smaller 
than r s (which is a reasonable assumption because all the 
analytical theories so far suggested that the PE persis- 
tence length scales as r s or r 2 ). Thus, the self energy of 
the chain can be written as E — Lri 2 h\(r s /lB}/D. At 
length scale smaller than r s , the polyelectrolyte behaves 
as a neutral chain under an uniform tension 



F = dE/dL = rj 2 \n(r./l B )/D 



(13) 



The average angle a bond vector makes with respect to 
the axis of the chain, therefore, is: 

fo exp(F/s cosO/ksT) cos#sin#d# 



coth 
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The charge density T) can be calculated as 

rj = rjo/ (cos 9) . (15) 



At weak screening r s 3> Ib, it can be estimated analyti- 
cally: 
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where the expansion terms of the order of l/ln 2 (r s /^) 
and higher were neglected. 

The more accurate numerical solution of Eq. ( |l3| ) , dl4| ) 
and (|l5|) for r\ is plotted in Fig. |^ by the thick solid line. 
One can see that the values rj(r^) calculated experimen- 
tally using OSF theory with growing N converge well 
to the theoretical curve for N = oo. Remarkably, the 
theoretical estimate for rj does not use any fitting pa- 
rameters. This, once again, strongly suggests the OSF 
theory is valid for flexible PE as well. 

The Bresler-Frenkel formula, Eq. ( ph, is also used to 
extract the persistence length in Ref. [L5| where the au- 
thors concluded that the dependence of l e on r s is sub- 
linear. The authors, however, used in Eq. (||) the bare 
contour length L, or in other words 77 ~ 770 , for the cal- 
culation of l e . As one can see from Fig. ||| this leads 
to 20-30% overestimation of the contour length of the 
effective chain where OSF theory is supposed to apply. 
To show that this overestimation is crucial, let us treat 
our data similarly to Ref. [l5| using rj = 770- We plot the 
resulting dependence of l e /r s on r s (similarly to Fig. 4 
of Ref. [15]) and compare it with our own results using 
corrected rj. The case N — 1024 is shown in Fig. |^. Ob- 
viously, the two results are different qualitatively. While 
the upper curve follows Eq. ([l2]) with slightly decreasing 
77, the lower curve shows sublinear growth of l e with r s 
(l e /r s is a decreasing function of r s ). This sublinear de- 
pendence observed in Ref. Il5| is clearly a manifestation 
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of their overestimation of the PE length L which should 
be used in Eq. @. 

Note that the true l e /r s curve should also eventually 
decrease to zero because l e saturates to the constant value 
r 1 2 L 2 /72Dk B T_ l when r s > L [See Eq. @]. But, accord- 
ing to Eq. 
where r s /L 



12]); this decay starts only at very large 
0.25. The deviation from l e oc r 2 at large 



r s seen in Fig. || is due to both the violation of the in- 
equality r s <C L and to the slight decrease of rj with r s . 



V. BOND ANGLE CORRELATION FUNCTION 

Another standard procedure used in literature is to cal- 
culate the persistence length of a polyelectrolyte as the 
typical decay length of the bond angle correlation func- 
tion (BACF) along the contour of the chain, assuming 
the later is exponential. 



f(\s'-s\) = (cos[Z(b s ,b s >)]) oc exp - 



It. 



(17) 



Here b s and b s / are the bond numbered s and s' respec- 
tively and Z(b s ,b s /) is the angle between them. The 
symbol (...) denotes the averaging over different chain 
conformations. To improve averaging, the pair s and s' 
are also allowed to move along the chain keeping \s' — s\ 
constant. 

We argue in this section that this method of determin- 
ing persistence length actually has a very limited range 
of applicability. At either small or large r s , the results of 
persistence length obtained from BACF are not reliable. 
In the range where this method is supposed to be ap- 
plicable, we show that the obtained L are close to those 



obtained in Sec. IV above. 

For small r s , excluded volume plays important role 
and, strictly speaking, it is not clear whether BACF is ex- 
ponential, and if yes, how one should eliminate excluded 
volume effect and extract l p from the decay length. Ac- 
cording to Ref. [l4|, the decay is not exponential in this 
regime. 

The procedure of determining the persistence length 
using BACF becomes unreliable at large r s as well. To 
elaborate this point, in Fig. ^a, we plot the logarithm 
of the bond angle correlation function f(x) along the PE 
contour length for a N = 512 and r s = 501b, typical 
values of N and r s where the excluded volume due to 
Coulomb interactions is small. There are three regions 
in this plot. In region A at very small distance along the 
PE contour length, monomers are within one electrostatic 
blobs from each other and the effects of Coulomb interac- 
tion are small. The bond angle correlation in this region 
decays over one bond length lg. At larger distance along 
the PE contour length, the region B, the decay is expo- 
nential and a constant decay length seems well-defined. 
Finally, at distance comparable to the chain's contour 
length, one again observes a fast drop of the BACF (re- 
gion C). This end effect is due to the fact that the stress 
at the end of the chain goes to zero and the end bonds 




FIG. 6: The logarithm of the bond correlation function f(x) 
as a function of the distance x (in units of is) along the chain 
for the case TV = 512, r s = 501b- There are three regions A, 
B and C. The dotted line, -0.47 - x/1083, is a linear fit of 
region B suggesting that the persistence length for this case 
is l p = 1083Z S . 



TABLE I: Comparison between Ibacf calculated using 
BACF method and r/le/rjo calculated in Sec. |Fy|. All lengths 
are measured in units of Is- 



N 


r s 


Ibacf 




2048 


100 


4590 


3682 




150 


9000 


7484 


1024 


80 


2535 


2180 




100 


3733 


3111 


512 


50 


1083 


809 



become uncorrelated. The persistence length of interest 
can be defined as the decay length in region B. 

Problem arises, however, at large enough r s when the 
region C (the end effect) becomes so large that region 
B is not well defined. In this case the obtained decay 
length underestimates the correct persistence length. As 
one can see from Fig. ^, region C can be quite large. It 
occupies 40% of the available range of x, even though the 
screening length is only 10% of the contour length in this 
case. 

There is an even more strict condition on how large r s 
is when the method of BACF loses its reliability. If l e is 
larger than L, the decrease of In f(x) in region B is less 
than unity. When this happens, an exponential decay is 
ambiguous. 

Because of all these limitations, in this section we use 
BACF to calculate l e only in the very limited range of 
r s where excluded volume is not important and l e is not 
much larger than L (the decrease in region B is greater 
than 0.1). The obtained Ibacf, which is measured along 
the chain contour, is compared to rjl e /r)o obtained using 
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the Bresler-Frankel formula in the previous subsection. 
(The factor 77/770 is needed because Ibacf is measured 
along the real PE contour while l e is measured along the 
renormalizcd PE contour.) The results are shown in the 
Table |. The two persistence lengths are within 20-25% 
of each other. This reasonably good agreement between 
two different methods shows that our calculations are 
consistent. It further strengthens the conclusion of two 
previous sections that OSF theory is correct in describing 
flexible polyelectrolytes. 

VI. CONCLUSION 

In this paper, we use extensive Monte Carlo simulation 
to study the dependence of the electrostatic persistence 
length of a polyelectrolyte on the screening radius of the 
solution. Not only did we simulate a much longer poly- 
electrolyte than those studied in previous simulations in 
order to show the scaling ranges, we also used a refined 
analyses which take into account local fluctuations to cal- 
culate the persistence length. These improvements result 
in a good support for OSF theory. They also help to 
explain why previous simulations failed to support OSF 
theory. 

In order to describe our numerical data we used a mod- 



ified OSF theory in the framework of ideas of KK. Linear 
charge density 77 was corrected to allow for short range 
fluctuations. In our case this is a relatively small correc- 
tion to 770 because we deal with a strongly charged PE. 
When one crosses over to sufficiently weakly charged PE 
linear charge density becomes strongly renormalized and 
matches KK expressions. We confirmed that corrections 
of 77 do not affect dependence of persistence length 
which was predicted by OSF for Iq <C r s <C L. In other 
words, we confirm KK idea that at large r s all effects 
of flexibility of PE are limited to a renormalization of 
77. At r s comparable to contour length L we found a 
good agreement of the numerical data with OSF formula 
modified for this case [Eq. (|l2|)] , which is derived in Ref. 
|l|. Again all effects of local flexibility are isolated in the 
small correction to the linear charge density 77. 
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